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Abstract: 

Effects of temperature dependent viscosity and thermal conductivity in a mixed convection boundary 
layer flow of a micropolar fluid towards a heated shrinking sheet in presence of magnetic field along 
with stagnation flow have been studied in the present work. The boundary layer equations are 
transformed into ordinary differential equations using similarity transformations. The effects of variable 
viscosity, variable thermal conductivity and the parameters involved in the study on the velocity, micro - 
rotation and temperature distribution profiles are investigated by solving the governing transformed 
ordinary differential equations with the help of Runge-Kutta 4 th order method with shooting technique. 
The numerical results are shown graphically and discussed in detail. 
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I. Introduction 

The problems of flow and heat transfer in the boundary layers of a continuous stretching/shrinking 
surface have attracted considerable attention of researchers due to their numerous applications in industrial 
manufacturing processes. Some of the applications are extraction of polymer sheets, paper production, hot 
rolling and glass-fiber production. Eringen [1] formulated the micropolar fluid theory as an extension of the 
Navier- Stokes model of classical hydrodynamics to facilitate the description of the fluids with complex 
molecules .The micropolar fluids are generally defined as isotropic, polar fluids in which deformation of 
molecules is neglected. Physically, a micropolar model can represent fluids whose molecules can rotate 
independently of the fluid stream flow and its local vortices. Micro polar fluids have important applications in 
colloidal fluids flow, blood flows, liquid crystals, lubricants and flow in capillaries, heat and mass exchangers 
etc. 

Stagnation point flows have also applications in blood flow problems, the aerodynamics extrusion of 
plastic sheets, boundary-layer along material handling conveyers, the cooling of an infinite metallic plate in a 
cooling bath, and textile and paper industries. Flows over the tips of rockets, aircrafts, submarines and oil ships 
are some instances of stagnation flow applications [2]. Hiemenz [3] started the study of stagnation flow problem 
and reduced the Navier-Stokes equations for the forced convection problem to an ordinary differential equation 
of third order by using similarity transformation. Chamkha [4] solved the problem of the laminar steady viscous 
flows near a stagnation point with heat generation/absorbing. The steady two dimensional point flow of a power 
law fluid over a stretched surface was studied by Mahapatra and Gupta [5]. The numerical solution of unsteady 
boundary-layer flow of an incompressible viscous fluid in the stagnation point region over a stretching sheet 
was presented by Nazar et.al [6] (using Keller box method). The problem of steady two dimensional laminar 
MHD mixed convection stagnation point flow with mass transfer over a heated permeable surface was examined 
by Abdelkhalek [7]. Two dimensional steady incompressible mixed convection non orthogonal stagnation flow 
towards a heated or cooled stretching vertical plate was considered by Yian et.al [8]. The solution of hydro 
magnetic steady laminar two dimensional stagnation flow of a viscous incompressible electrically conducting 
fluid of variable thermal conductivity over a stretching sheet was obtained by Sharma and Shing [9] using 
shooting method. Effect of viscous dissipation on heat transfer in a non-Newtonian liquid film over an unsteady 
stretching sheet was investigated by Chen [10]. Heat transfer over a stretching surface with uniform or variable 
heat flux in micropolar fluids was studied by Ishak et. Al [1 1]. 
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Mahanta et .at. [12] discussed the effects of radiation of flow of second grade fluid over a stretching 
sheet through porous medium with temperature dependent viscosity and thermal conductivity. A numerical 
study of steady incompressible micro polar fluid in a two dimensional stagnation point flow towards a stretching 
sheet was investigated by Nazar et.al [13]. Chang [14] presented the solution of flow and heat transfer 
characteristic of mixed convection in a micropolar fluid along a vertical flat plate with conduction effects. Ishak 
et.al [15] studied steady stagnation flow towards a vertical surface immersed in a micropolar fluid. Rashidi et.al 
[16] investigated the mixed convection boundary-layer flow of a micropolar fluid towards a heated shrinking 
sheet. 

In the above literatures, in most of the studies, the viscosity and the thermal conductivity of the ambient 
fluid were assumed to be constant. When the effects of temperature dependent viscosity and thermal 
conductivity are taken into account, the flow characteristics are significantly changed compared to the constant 
property case. In this paper, an attempt has been made in this study to find the effects of temperature dependent 
viscosity and thermal conductivity on a mixed convection two dimensional stagnation point flow and heat 
transfer of a steady viscous incompressible micropolar fluid towards a heated shrinking sheet in presence of 
magnetic field with viscous dissipation. Viscosity and thermal conductivity are assumed to be inverse linear 
functions of temperature. 



II. Formulation of the Problem: 

We consider two dimensional stagnation point flow of a micro polar fluid impinging normally on a 

heated shrinking sheet at a fixed flat plat coinciding with the plane y = 0 . The flow is assumed to be laminar, 

steady, viscous and incompressible and except the fluid viscosity and thermal conductivity all the fluid 
properties are assumed to be constant. Also a magnetic field of constant intensity is assumed to be applied 
normal to the surface and the electrical conductivity of the fluid is assumed to be small so that the induced 
magnetic field can be neglected in comparison to the applied magnetic field. Under these assumptions the 
governing equations of the problem are as below: 
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Where u , v are velocity components in the directions of x and y along and perpendicular to the surface 
respectively . N is the component of micro -rotation vector normal to the xy -plane, p is the density , k is 
the vortex viscosity, g is the acceleration due to gravity, P is the coefficient of thermal expansion , y is the 
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spin gradient viscosity , j is the micro-inertia density , p is the pressure of the fluid , T is the temperature , 
v and 2 are the viscosity and thermal conductivity respectively which are the functions of v and y . c is 

the specific heat capacity at constant pressure of the fluid , B is the magnetic intensity and cr is the electrical 
conductivity . The term ± g j3 (T - T ) of equation (2) indicates the buoyancy force, where “+” sign refers 
buoyancy assisting and “-’’sign corresponds to the buoyancy opposing the flow regions. 

The boundary conditions for the problem are: 



u(x, 0) = bx, v(x,0) = 0, N (x, 0) = 0, 



T(x, 0) = T 



> 



(5) 



u ( x , oo ) = U = ax, 



N (x, oo ) = 0 



T (x,oo) = T 

\ 7 s r, 



Where b < 0 for the shrinking sheet, T w is temperature on the surface, T ^ is temperature of the fluid at infinity 
and U is the free stream velocity of the fluid. 

Following Lai and Kulacki [17] let us assume that, 
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where is the viscosity at infinity, A and T r are constants , T r is transformed reference temperature , 5 is a 

constant based on thermal property of the fluid and A < 0 for gas , A > 0 for liquid. Similarly, s and T are 
constants and their values depend on the reference state and thermal properties of the fluid i.e., 

To solve equations (1) - (4) subject to the boundary conditions given in equation (5) we use the following 
similarity transformations, 



n = ( — ) 2 y , p(x,co) = P - — — (x 2 + y 2 ) , u(x,y) — axf ( 77 ), 
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Where 77 is the similarity parameter, p is the stagnation pressure. Also from equations ( 6 ) and (7), we have, 
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Equation of continuity in equation (1) is identically satisfied using equation ( 6 ) and therefore the velocity field is 
compatible with continuity equation and represents the possible fluid motion. Using equations (7) and ( 8 ) in 
equations (2)-(4) we get the following differential equations: 
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6 r and 6 are the dimensionless parameters characterizing the influence of viscosity and thermal conductivity 
respectively, can be written as: 



6 




6 








y 



( 12 ) 



The transformed boundary conditions are, 



f(ri) = 0 , 
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Here the case a = 0 stands for Hiemenz flow towards a solid plate and the case oc > 0 is for the stagnation 
point over a stretching sheet. In our case of stagnation flow towards a shrinking sheet, hence we take oc < 0. 

The two important physical quantities of our interest in the problem are skin friction 
coefficient (cy) and Nusselt number (Nu) which are defined as, 
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III. Results and Discussion: 

The system of differential equations (9 - 11) together with the boundary conditions (13) is solved 
numerically by using the fourth order of Runge-Kutta integration accompanied with the shooting iteration 
scheme. The purpose of this study is to bring out the effects of the variable viscosity and variable thermal 
conductivity on the governing flow with the combinations of the other flow parameters. 

The numerical computations have been carried out for various values of magnetic parameter (M), Prandtl 
number (Pr ), Eckert number ( E ), buoyancy parameter ( R ), coupling constant parameter ( K ), the variable 

viscosity parameter (0 r ), variable thermal conductivity parameter (0 ) and micro inertia density parameter 

( A ) . In order to illustrate the results graphically, the numerical values of dimensionless velocity distribution 

/ ' ( rj ) , dimensionless micro-rotation distribution g ( 77 ) and temperature distribution 6 ( rj ) with variation of 
different parameters have been plotted in Figures 1-12. 




Figure 1. Variation of / ' (77 ) for different values of 6 
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Figure 2. Variation of / ' ( tj ) for different values of K 




Figure 3. Variation of ffrj ) for different values of M 




Figure 4. Variation of ffrj ) for different values of R 
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Figure 5. Variation of g (tj ) for different values of 6 




Figure 6. Variation of g (?j ) for different values of K 




Figure 7. Variation of g (77 ) for different values of M 




Figure 8. Variation of g (77 ) for different values of A 
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Figure 9. Variation of g (77 ) for different values of R 




Figure 10. Variation of 6 (77 ) for different values of 0 c 
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Figure 11. Variation of 6 (77 ) for different values of M 




Figure 12. Variation of 6 (77 ) for different values of E 
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The Figures 1-4 represent the velocity distribution with the variation of viscosity parameter 6 r coupling 

constant parameter K , magnetic parameter M , Buoyancy parameter R respectively. It is seen that velocity 
increases with the increasing values of 9 r and R where as it decreases with the increasing values of M and 

K . It has been observed that with the variations of thermal conductivity parameter# , micro-inertia density 

parameter A , Eckert number E and Prandtl number Pr the variation of velocity is not significant. 

Figures 5-9 represents the variation in micro-rotation distribution with the variation of viscosity parameter 
6 r , coupling constant parameter K , magnetic parameter M , micro-inertia density parameter A and 

Buoyancy parameter R respectively. It is observed that micro -rotation decreases with the increasing values of 
0 r and R respectively. It is evident from the Figures 6, 7 and 8, the micro-rotation distribution increases with 

the increasing values of K , M and A respectively. In this case also variation of micro -rotation is not so much 
significant with the variations of thermal conductivity parameter 9 , Eckert number E and Prandtl number P r . 
The Figures (10-12) represents variation of temperature distribution with the variation of thermal conductivity 
parameter 9 , Magnetic parameter M and Eckert number E . It is observed that temperature decreases as 9 

increases whereas it increases with the increasing values of M and E . In several practical applications, the 
surface characteristics such as friction factor and Nusselt number play important roles and hence, the missing 
values of / "(0) , g'(0) and #'( 0) for various values of 9 r ,9 c9 M ,Pr, K and E have been derived in 
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From the Table- 1, it is observed that with the increasing values of 0 r and E c , the values of / "( 0) and g'(0) 

are increasing and the values of O’ ( 0) are decreasing for the increasing values of 6 r when all other parameters 

are fixed. It is seen that with the increasing values of 0 c andPr , the values of / "(0) , g '(0) and 0\ 0) are 

decreasing. It is also indicated that when the magnetic parameter M and coupling constant parameter K 
increases, the values of 6 '( 0 ) increases but the reversal trend is observed for /"( 0) and g'(0) 

IV Conclusions 

In this study , the effects of temperature dependent viscosity and thermal conductivity in a mixed 
convection boundary layer flow of a micropolar fluid towards a heated shrinking sheet in presence of magnetic 
field is examined. Numerical solutions are presented for the fluid flow and heat transfer characteristics for 
different values of parameters involved in the problem. The effects of temperature dependent viscosity and 
thermal conductivity on velocity, micro -rotation and temperature distribution are quite significant. Thus, the 
present study will serve as a scientific tool for understanding more complex flow problems concerning with the 
various physical parameters. 
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